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Abstract. We study how Type la/II supernova rates at 
various redshifts can be used to constrain the cosmic star 
formation rate (CSFR) history. We use a spectrophoto- 
metric model of galaxy evolution which provides a self- 
consistent description of global properties of local galax- 
ies. We have chosen two CSFR histories, with and without 
dust extinction, both normalized in order to reproduce the 
observed comoving luminosity densities at different wave- 
lengths. We show that as expected, the SNII rates pro- 
vide a direct measurement of the instantaneous CSFR, 
whereas the SNIa rate depends on the CSFR history dur- 
ing the past few Gyrs. Therefore the SNIa rate detected 
up to z~ 1 is able to provide new insight on the CSFR 
at higher redshift. The comparison with the data are in 
good agreement with both local and high-z measurement 
and provides an additional evidence of the validity of the 
CSFR up to z = 1, but suggests values at z > 1 higher 
than those of Madau et al. (1996) and consistent with the 
presence of extinction. We also show that Type la super- 
nova observations at Z'^ 1 can already put constraints on 
the still puzzling nature of Type la progenitors. 



1. Introduction 

Supernova (SN) rates at high redshift are very important 
not only because of their ability to allow the direct de- 
termination of cosmological parameters such as Ho, f^o 
and r^A, but also in understanding nucleosynthesis rates, 
galaxy evolution and star formation. More specifically, 
the evolution of Type la SN rates with redshift can al- 
low one to probe the past history of star formation in 
the universe, whereas their Type II SN counterparts di- 
rectly probe the instantaneous star formation because of 
the short-lived massive stars which explode. Furthermore, 
observations of high-z Type la supernovae would help us 
to better understand the nature of the supernova progen- 
itors (Ruiz-Lapuente and Canal 1998). On the observa- 
tional side, there has been a renewed interest in the search 
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of SNe (Cappellaro et al. 1997 hereafter C97). An increas- 
ing number of Type la SNe, which are the most homo- 
geneous and most luminous SNe, are detected at redshift 
z ^ 1 (Perlmutter 1997, Tonry et al. 1997) and recently 
the Supernova Cosmology Project has reported the first 
measurement of the Type la SN rate at z^ 0.4 (Pain et 
al. 1997; hereafter P97). Thanks to recent deep redshift 
surveys such as the Canada-France Survey (CFRS) and 
the Hubble Deep Field (HDF), it has been possible to 
study galaxy populations as a "whole" in order to infer 
the global history of star formation in the universe. More 
specifically, Madau et al. (1996; hereafter M96) have inves- 
tigated the cosmic star and metal formation in the early 
universe by combining ground-based data with the HDF 
observations. In this letter we use the chemical and photo- 
metric population synthesis code of Sadat and Guiderdoni 
(1998; hereafter SG98) to make predictions on the cosmic 
evolution of SN rates with redshift based on the CSFR de- 
rived from observations. We expect that such studies can 
provide additional constraints on the CSFR history and 
metallicity evolution and will allow one to interpret SN 
rate measurements at high redshift which are now becom- 
ing available. We hereafter adopt Hq = 50 kms~^Mpc^^, 
57o = 1 and r^A — 0. 

2. Type la and II Supernovae 

The model we use is a spectrophotometric model based 
on stellar libraries from Kurucz (1992) supplemented by 
Bessell et al. (1989, 1991) for M Giants and Brett (1995) 
for M dwarfs. The stars are followed along stellar tracks 
computed by the Geneva team (Charbonnel et al. 1996). 
The Initial Mass Function we choose is a "standard" 
IMF defined by: (/)(to) oc m~^, with index a;=0.25 for 
0.1 < m < IMq, 1.35 for 1 < m < 2Mq and x =1.7 
for 2 < m < 12OM0. A more detailed description of the 
model is presented in a forthcoming paper (SG98). Note 
that, as a first step, we here restrict ourselves to models 
with a single metallicity (Z=Zo). The chemical evolution 
will be studied in SG98. The progenitors of Type II,Ib,c 
SNe are stars with masses M > 8Mq. The rates are di- 
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Fig. 1. Evolution with redshift of the luminosity density at rest-frame wavelengths. The data points with error bars are taken 
from Lilly et al. (1996) {filled squares), and Connolly et al. (1997) {empty hexagons). Open triangles and filled triangles respec- 
tively correspond to 1600 A HDF measurements uncorrected (Madau 1997) and corrected (Pettini et al. 1997) for extinction. 
The empty square is from Treyer et al. (1997). The two curves correspond to the two different parameterizations Ml {solid line) 
and M2 {dashed line) of the CSFR. 



rectly computed from the IMF. The nature of the Type 
la SNe progenitors is still a matter of debate. Indeed, it 
is not yet clear which of the two main competing models 
- double degenerate (DD, Iben & Tutukov 1984) or sin- 
gle degenerate (SD, Whelan & Iben 1973) - applies for 
the SNIa precursors. Here we assume the SD model where 
Type la SNe are produced by C-ignition and total disrup- 
tion of a cold degenerate WD when this latter exceeds the 
Chandrasekhar mass after mass transfer in a close binary 
system. We estimate the rate of Type la SNe according to 
the formalism by Ferrini et al. (1992; see also Greggio & 
Renzini 1983). In our SD modelling, the free parameters 
'm-Binf and ao) are respectively the lower limit of the total 
mass of the binary systems which can produce Type la 
SN, and the fraction of the total mass of stars which be- 
long to these systems. They have been adjusted in order to 
reproduce the main properties of the solar neighbourhood. 
We find that the couple [rnBm/,ao] — [3,0.05] correctly 
reproduces the data. In order to assess the uncertainties 
in the nature of the progenitors, we have also used the 
empirical parameterization of the rates of SN la as de- 
rived by Ciotti et al. (1991; hereafter C91) Rja on ^sArt^". 
9sN is the normalizing parameter. The Rja evolution is 
controlled by the free parameter s. In order to account 
for the Fe content in clusters of galaxies, C91 concluded 
that s > 1.4. Note that our standard modelling roughly 
corresponds to s '-~^ 1.6. The other free parameter is the 
rise time which is fixed to ti5.o=0.05 in units of 15 Gyr 
(Renzini et al. 1993). 



3. The Cosmic Star Formation Rate 



By combining different photometric surveys such as the 
CFRS, up to ;z = 1 (Lilly et al. 1995) and the HDF up 
to z = 4 (M96), it has been possible to derive a picture 
of the star formation rate history of the whole universe 
(M96, Madau et al. 1997; M97). To compute the evolu- 
tion of the cosmic supernova rates per comoving volume 
(CSNR), we make use of the CSFR as derived in M97 
(model Ml). We also introduce another CSFR law to ac- 
count for a possible dust extinction correction (model M2). 
Ml and M2 have the same shape at low z in agreement 
with the CSFR derived from observations, but at higher z, 
we have allowed the CSFR to be higher in model M2 than 
in Ml in agreement with Pettini et al. (1997). These time- 
dependent CSFRs are then introduced as an input in our 
code. The amplitude is normalized in order to match the 
luminosity density evolution in the UV-continuum with 
the following constants of proportionality between the ob- 
served Luv and our derived SFR: (7.5 x 10^^,6.5 x 10^^) 
in units of W Hz-i (Mq yr"!)"! at (1500 A, 2800 A). Be- 
fore computing the SN rates, we first have checked whether 
our simple model with a "standard" IMF and the param- 
eterization of the CSFR is able to reproduce the observed 
volume averaged density luminosities p^, at different A. 
Figure 1 shows that the agreement at longer A is satisfy- 
ing. The 4400 and 10000 A luminosity densities seem to 
hint at model M2 in agreement with the IR/submm back- 
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Fig. 2. Evolution of SN rest-frame rates per comoving volume 
for model Ml {lower curves) and M2 {upper curves): Ru/ib.c 
{solid lines); Ria {dashed lines). CSFRs are shown as dotted 
hues. Data are derived from the UV fluxes with the constants 
of proportionality mentionned in section 3. 
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Fig. 3. Predicted evolution of the rest-frame Type la SN rate 
per comoving volume. The curves labelled with s correspond to 
Ria as derived assuming the C91 parameterization normalized 
in order to reproduce the present day value for the same CSFR. 
The horizontal bar is the observational rcdshift range. 



ground which also suggests the presence of extinction at 
high z (Guiderdoni et al. 1997). 

4. Cosmic Evolution of SN Rates 

In this section, we use the CSFR to derive the CSNR for 
both II/Ib,c and la Types. To our knowledge, this is the 
first time that the evolution of SN rates with redshift is 
predicted from a self-consistent spectro-photometric mod- 
elling of galaxy evolution at high z and independently of 
the details of individual galaxy evolution. As already men- 
tionned, the direct measurement of SNe can be used as an 
independent test for the cosmic star and metal formation 
in the universe. In figure 2, we have plotted the predicted 
evolution of the CSNR per unit of comoving volume with 
redshift. The Type II/Ib,c rate shows the same shape as 
the instantaneous CSFR, that is the rise, peak and drop 
from high redshift to the present time. This means that 
the SNII/Ib,c rate can be used as an independent tracer of 
the star formation rate. The Type la rate Ria has a differ- 
ent shape from the CSFR while coincidcntly it has nearly 
the same behaviour as the B-band luminosity. The most 
important point is the time delay we observe between the 
stellar (binary in this case) birth and the explosion time. 
The occurrence of the SNIa rate peak is shifted by a few 
Gyrs relatively to the CSFR peak. However, at z > 0.9, 
Type la SN can be used as a probe of the past history of 
the CSFR. Furthermore, as can be seen from figure 3, Ria 
has very distinct shape (different normalization and dif- 
ferent time of peak occurrence) depending on the adopted 
CSFR shape, and this difference is higher at higher z. 
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Fig. 4. Evolution of SNII {upper curves) and SNIa {lower 
curves) rates for models Ml{solid lines) and M2 {dotted lines) 
compared to the observations. The filled circle and open circle 
respectively correspond to SNII/Ib,c and SNIa rates from C97, 
and the filled triangle is from P97. 



Therefore, measurements of SNIa at z^ 1 would be able 
to discriminate between the models. 

5. Comparison with observations 

Supernova rate has been measured in nearby galaxies by 
several authors (Cappellaro & Turatto 1988, Evans et al. 
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1989, C97). At higher redshift, Pain and his collaborators 
have recently reported the first high-z (z ^ 0.4) rest-frame 
Type la SN rate with a value of -- 0.82 h^ SNu. In order 
to compare the model to the observations, we have con- 
verted the observed rate in SNu (SNe/100 yr/ 10^" Lbq) 
into a rate in yr~^ using blue luminosities as computed 
by our code. From figure 3, we can see that although the 
available measurement of the SNIa rate does not allow one 
to discriminate between models Ml and M2, the redshift 
is still low and the statitics is poor. Measurement at z ^ 
1 with the same statistics, would begin to discriminate 
between models, independently of the adopted model for 
Type la progenitors. By investigating different models for 
SNIa rates using C91 empirical parameterization (see fig- 
ure 3 caption), wc have found that a comparison between 
these model predictions and observations at low and in- 
termediate z shows that the best match is obtained for 
models with a rapidly dechning Type la rate (s > 1.4) 
consistent with the constraints on s from chemical evo- 
lution. In figure 4 we have plotted the predicted SNIa/II 
rates in SNu. As can be seen, the SNIa rate is nearly con- 
stant. This is expected since both the RsNia per unit 
of comoving volume and the B-band luminosity have the 
same shape. Therefore the peak observed in RsNia has 
been erased, leading to a constant ratio. However, this is 
not true for the Type II SN rate which turns out to be 
a monotonic function increasing with redshift. This is a 
consequence of the fact that the stellar population which 
emits the bulk of the B-band luminosity is older than the 
SNII progenitors. Note that the SN rates in SNu do note 
change drastically with the adopted CSFR prescription. 
This means that SN rates in SNu loose a large part of 
the information on the absolute value of the CSFR. Our 
model reproduces fairly well the observed local Type la 
rate and our prediction for high-z rates is in good agree- 
ment (within the error bars) with the observed Type la 
rate of P97. At z ^ 1 we expect that the SNIa rate in SNu 
does not change from its local value while the SNII rate 
is expected to be 3 times the local value. This leads to 
the decrease of the Type la/II SN ratio from z=0 to z=l 
(see also Yungelson & Livio 1997). However, the Type II 
SN rate is not observed yet and this prediction could be 
confirmed with future instruments such as the NGST. 



6. Conclusion 

We have computed the CSNR from the CSFR using a 
self-consistent modelling. We found that: 

— The adopted standard IMF allows one to fit the ob- 
served local and high-z colours of the universe. 

— Type II SN rates are very good tracers of instantaneous 
star formation and could be used to set constraints on 
the CSFR at 2 < z < 4 with the new generation of 
telescopes such as the NGST. 



— Type II SN rates Rsnii (expressed in SNu) are 
an increasing function with redshift while the 
RsNia/ Rsnii ratio is decreasing. 

— Our model predictions for SN rates are in a good agree- 
ment with observations locally as well as at higher z, 
and consistent with the current limits on the CSFR. 

— At redshift z ^ 1 , a comparison between the predicted 
and the observed Type la rates would allow us to put 
stong constraints on models for SN la progenitors. In 
the frame of C91 models, wc found that the best match 
is obtained for models with a steep evolution of SNIa 
rate (i.e. s > 1.4) consistent with the chemical evolu- 
tion of Fe in clusters of galaxies. 

— Intermediate and high-z Type la SN rates can be used 
as an independent test to constrain the CSFR. The 
current data at z ~ 0.4 are not sufficient to disentan- 
gle between the models but observations of Type la 
supernovae at z ^ 1 would be critical for understand- 
ing the CSFR and would allow one to assess whether 
the z > 1 CSFR is higher than directly observed from 
the UV. 
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